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CHAPTER THREE 
TWO DIMENSIONAL WAVES 



1 Reflection and tranmission of sound at an inter- 
face 

Reference : L. M. Brekhovskikh and O. A. Godin: Acoustics of Layered Media I. 
Springer. §.2.2. 

The governing equation for sound in a honmogeneous fluid is given by (7.31) and 
(7.32) in Chapter One. In term of the the veloctiy potential defined by 

u = V0 (1.1) 

it is 

^It^VV (1.2, 
where c denotes the sound speed. Recall that the fluid pressure 

p = -pd(p/dt (1.3) 

also satisfies the same equation. 

1.1 Plane wave in Infinite space 

Let us first consider a plane sinusoidal wave in three dimensional space 

0(x, t) = (/)oe^^'*-^-'"*) = 0oe^('="-^-'"*) (1.4) 



Here the phase function is 

^(x, t) —\s.-'x. — ut 



(1.5) 
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The equation of constant phase 0{:x.,t) — 9o describes a moving surface. The wave 
number vector k = /cn is defined to be 

k^kn^Ve (1.6) 

hence is orthogonal to the surface of constant phase, and represens the direction of wave 
propagation. The frequency is defined to be 

Is (1.4) a solution? Let us check (1.2). 

\ox ay oz J 



VV = V ■ V0 = ik ■ ikcj) = -k\ 



5V 



Hence (1.2) is satisfied if 

u = kc. (1.8) 
Sound in an infinite space is non-dispersive. 



1.2 Two-dimensional reflection from a plane interface 

REfering to figure ??, let us cnsider two semi-infinite fluids separated by the plane 
interface along z — 0. The lower fluid is distinguished from the upper fluid by the 
subscript " 1" . The densities and sound speeds in the upper and lower fluids are p, c and 
Pi, Ci respectively. Let a plane incident wave arive from 2; > at the incident angle of 
9 with respect to the z axis, the sound pressure and the velocity potential are 

Pi — Pq exp[ik{x sinO — z cos 9)] (1.9) 

The velocity potential is 

iP 

(f)i — exp[ik{x sm9 — z cos 9] (1-10) 

ujp 

The indient wave number vector is 



k^ = {ki,ki) = k{sm9,-cos9) 



(1.11) 
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Figure by MIT OCW. 

Figure 1: Plane wave incident towards the interface of two fluids. 

The motion is confined in the x, z plane. 

On the same (incidence) side of the interface we have the reflected wave 

Pr = -Rexp[iA;(xsin6' + 2; cos 6*)] 

where R denotes the reflection coefficient. The wavenumber vector is 

k^' = {kl,kl) = A:(sin^,cos^) 

The total pressure and potential are 

p = Pq {exp[ik{xsin6 — 2;cos^)] + i?exp[iA;(xsin^ + ^cos^)]} 
iP 

= {exp[2/i;(xsin6' — 2;cos6')] + Rexp[ik{x sin 9 + 2; cos 6')]} 

pu 

In the lower medium z < the transmitted wave has the pressure 

Pi = TPoexp[iA;i(xsin6'i — zcos6'i)] 
where T is the transmission coefficient, and the potential 



-Texp[zA;i(xsin^i — 2;cos^^i)] 



:i.l2) 

:i.i3) 
:i.i4) 

:i.l5) 
:i.l6) 

:i.l7) 



Along the interface z = we require the continutiy of pressure and normal velocity, 



I.e., 



p = Pi, z = 



:i.i8) 



3.1. Reflection of sound by an interface 



4 



and 

w^wi, z = 0, (1.19) 

Applying (1.18), we get 

|gifcxsme ^ ^gifcxsinej _ j. p^^ik^x sin 9, ^ - OO < X < OO. 

Clearly we must have 

ksine ^ kisinOi (1.20) 

or, 

sin 9 sin 9i 



C Ci 

This is just Snell's law. With (1.20), we must have 



Applying (1.19), we have 



— r-A;cos^e**^^'"^ + i?A;cos^e*'=^'"n = — l-h cos OiTe''"' 
pcu piUJ 



which implies 



Eqs (1.22) and (1.23) can be solved to give 



p\k cos 9 + pki cos 9 



1 



Alternatively, we have 



pc cos 9 1 + piCi cos 9 
piCi cos 9 — pc cos ^1 
PiCi cos 9 + pc cos ^1 



Let 



Pi c 



(1.21) 



l + i? = T (1.22) 



„ pki cos ^1 _ , , 

l-i?=^-^ (1.23) 



T = '^Pikcos9 

p/ci cos 9i + pi A; cos 6' 

P _ pikcos9 - pki cos6'i , , 



_ 2piCiCos^ 

T = (1.26) 



piCi cos 61 -pc cos 6li 



:i.28) 
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where the ratio of sound speeds n is called the index of refraction. We get after using 
Snell's law that 



mcos^-ncos^i mcos9-nJl 



R = = y ^ Q 29) 

mcos^ + ncos^i mcos^ + n./T^^ 



The transmission coefficient is 



^ 2m cos 6 

T=l + R= (^-2°) 

mcos6' + nA/l - 



We now examine the physics. 

1. If n = c/ci > 1, the wave enter from a faster to a slower medium, then 9 > 9i and 

sin^/?i < 1 always. R is real. In particular, for normal incidence 9 = 9i = 

m — n , ^ 

R= 1.31 

m + n 

is real. If m > n > 1, i.e., pi/p > c/ci > 1, then R is positive. If n > m > 1 then 
R is negative. As 9 increases, 9i also increases. R decreases until at 6* = 7r/2, so 
that 

R=-- = -l (1.32) 

n 

Hence R lies on a segment of the real axis as shown in Figur l.a. or Figure l.b. 
Note that there is no reflection i? = at the special incidence angle 9 = 9b, called 
the Brewster angle: 

sm9B = ^== (1.33) 

2. If however n < 1, i.e., the wave enters from a slower medium to a faster medium, 
then 9i > 9. for sufficiently small 9, R is real. When 9 increases to a critical value 
6, defined by 

sin5 = n (1.34) 

sin^i = 1 so that 9i becomes 7r/2, and 

m cos S 



R 



m cos 6 



Below this critical angle (0 < < 5), R is real. In particular, when 9 = 0, (1.31) 
applies. See Figure 2.c for m > n and in 2.d. for m < n. 
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Fig. 2.2a-d. Reflection coefficient in 
the complex plane for various rela- 
tions between the media parameters: 
3 = Sia) i<c/ci<ei/e; 

(b) 1 < c/ci, gi/g < c/ci; 

(c) c/ci < 1, c/ci < ei/g; 

(d) Qi/g < c/ci < 1 



Figure 2: Complex reflection coefficient. From Brekhovskikh and Godin §.2.2. 
When 9 > 5, the square roots above become imaginary. We must then take 

cos^i = 



, sin^^ . sin'^e , 
1 IT- = i\l — ^ 1 



1.35) 



The reflection coefficient is now complex 

mcos^-m^^ ^ 



mcose + inx/'^ - 1 



;i.36) 



with \R\ — 1, implying complete reflection. As a check the transmitted wave is 
now given by 



Pt^T exp 



ki I ixsin^i + z^^ sin^ — 1 



(1.37) 



so the amplitude attenuates exponentially in z as z — > — oo. Thus the wave train 
cannot penetrate much below the interface. The dependence of it! on various 
parameters is best displayed in the complex plane R — ^R + i^R. It is clear 
from (1.36 ) that QR < so that R falls on the half circle in the lower half of the 
complex plane as shown in Figures 2.c and 2.d. 
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2 Equations for elastic waves 

Refs: 

Graff: Wave Motion in Elastic Solids 

Aki & Richards Quantitative Seismology, V. 1. 

Achenbach. Wave Propagation in Elastic Solids 

Let the displacement vector at a point Xj and time t be denoted by Ui{xj,t), then 
Newton's law applied to an material element of unit volume reads 

d'^Ui dTij 

Pm^ = d^ ^2.1) 

where Tij is the stress tensor. We have neglected body force such as gravity. For a 
homogeneous and isotropic elastic solid, we have the following relation between stress 
and strain 

Tij = \ekkSij + 2(xeij (2.2) 
where A and /i are Lame constants and 

e.-UI:^^'^ (2.3) 



2 \dxj dx 

is the strain tensor. Eq. (2.2) can be inverted to give 



where 



is Young's modulus and 



e,^^r,-^r,kS, (2.4) 



E = ^^^^±J^ (2.5) 
A + /i 

A 



2(A + /i)' 

Poisson's ratio. 

Substituting (2.2) and (2.3) into (2.1) we get 

dxj dxj ^dxj \dxj dxi 
= \- 1^^ — ^ ^ A* 



(2.6) 



dxi dxjdxj dxidxj 
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In vector form (2.1) becomes 



= (A + /x) V(V • u) + /xV^u (2.7) 



Taking the divergence of (2.1) and denoting the dilatation by 



dui du2 dus ,^ 
ox I 0x2 0x3 



we get the equation governing the dilatation alone 



(A + /x) V • VA + /xV'A = (A + 2/x) V'A (2.9) 



or, 



where 



ciV^A (2.10) 



Thus the dilatation propagates as a wave at the speed c^. To be explained shortly, this 
is a longitudinal waves, hence the subscript L. On the other hand, taking the curl of 
(2.7) and denoting by cu the rotation vector: 

cJ = Vxu (2.12) 

we then get the governing equation for the rotation alone 

= 4V^a; (2.13) 



where 

CT^.f^ (2.14) 

V p 

Thus the rotation propagates as a wave at the slower speed ct- The subscript T indicates 
that this is a transverse wave, to be shown later. 

The ratio of two wave speeds is 



Ct V A* 

Since 



'^ + ''>l. (2.15) 
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it follows that the speed ratio depends only on Poisson's ratio 



Cl 
Ct 



2-2v 
l-2u 



(2.17) 



There is a general theorem due to Helmholtz that any vector can be expressed as 
the sum of an irrotational vector and a solenoidal vector i.e., 



subject to the constraint that 



u = V0 + V X H 



V-H = 



(2.18) 



(2.19) 



The scalar and the vector H are called the displacement potentials. Substituting this 
into (2.7), we get 

p— [V0 + V X H] = //V'[V(/) + V X H] + (A + //)VV • [V(/) + V x H] 



Since V • V0 = VV, and V • V x H = we get 



(A + 2/x) VV - P 



92 



+ V X 



(2.20) 



Clearly the above equation is satisied if 



(A + 2At) - p 



and 



AfV'H - p 







(2.21) 



(2.22) 



Although the governing equations are simplified, the two potentials are usually coupled 
by boundary conditions, unless the physical domain is infinite. 
A typical seismic record is shown in Figure 2. 



3 Free waves in infinite space 

The dilatational wave equation admits a plane sinusoidal wave solution: 



(3.1) 
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Waves in Elastic Solids 




ERM 

A typical seismogram during an eartliqual<e, sliowing 
tlie arrival of the P(primary), S(secondary) and R(Rayleigh) 
waves. Note that the Rayleigh waves have the largest amplitude. 



Figure by MIT OCW. 

Figure 3: A typical seismic record. P: Longitudinal wave; S: Transverse wave; R : 
Rayleigh surface wave. From Billingham & King. Wave Motion, Cambridge U Press. 
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Here the phase function is 

9{x,t)^k{n-x-CLt) (3.2) 
which describes a moving surface. The wave number vector k = /en is defined to be 

k = A;n = (3.3) 

hence is orthogonal to the surface of constant phase, and represents the direction of 
wave propagation. The frequency is 

dd 

uJ = kcT^-^ (3.4) 

A general solution is 

(l) = (pin ■ X — cit) (3.5) 
Similarly the the following sinusoidal wave is a solution to the shear wave equation; 

H = Hoe**^^""'-'^^*) (3.6) 

A general solution is 

H = H(n • X - CTt) (3.7) 



Note: 

We can also write (3.5) and (3.9) as 



cl 

and 



, / n • X, , , 

(j>^(j>{t ) (3.8) 



H = H(i - — ) (3.9) 

Ct 

where 

n n , , 

Sl = — , St = — 3.10 

Cl Ct 

are called the slowness vectors of longitudinal and transverse waves respectively. 



In a dilatational wave the displacement vector is parallel to the wave number vector: 

Ui = V0 = 0'n (3.11) 
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from (3.5), where 0' is the ordinary derivative of with repect to its argument. Hence 
the dilatational wave is a longitudinal (compression) wave. On the other hand in a 
rotational wave the displacement vector is perpendicular to the wave number vector, 
„ fdHz dHy\ fdH^ OH A fdHy OH. 



dy dz J \ dz dx J \ dx dy 

= n X H' (3.12) 
from (3.7). Hence a rotational wave is a transverse (shear) wave. 

4 Elastic waves in a plane 

Refs. Graff, Achenbach, 

Aki and Richards : Quantitative Seismology, v.l 

Let us examine waves propagating in the vertical plane oix,y. All physical quantities 
are assumed to be uniform in the direction of z, hence d/dz — 0, then 

^ ^dl^dHz ^ ^dl_dH^ ^ ^_dH^^dHy 
^ dx dy ^ ^ dy dx ^ ^ dy dx 

and 
where 

^ + ^ = 1^ (4 15) 

dx'^ dy"^ c\ dt"^ ' 

dx^ ^ dy^ 4 dt^ ' ^ ^'^'^ ^ ^ 

Note that Uz is also governed by (4.16). 

Note that the in-plane displacements Ux,Uy depend only on 4> and i?^, and not on 
Hx,Hy. Out-of-plane motion Uz depends on Hx,Hy but not on Hz- Hence the in- 
plane displacement components independent of the out-of-plane component 
Uz- The in-plane displacements (u^jUy) are associated with dilatation and in-plane 
shear, represented respectively by and Hz, which will be refered to as the P wave and 
the SV wave. The out-of-plane displacement Uz is associated with and Hy, and will 
be refered to as the SH wave. 
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Prom Hooke's law the stress components can be written 

f du^ duy\ , ^ <9m^ /X , o ^ /''^"^ , o ^^y 

T.. = A — + ^ +2/X— = (A + 2/x) — + ^ -2/.^ 



dx dy J dx \ dx dy J dy 



dx dy J dy \ dx dy J dx 

^^y-^' l^ + ^)=^'^7^-^ + ^] (4-20) 



c?x dy J \ dxdy dx"^ dy"^ 
du, ( d^H d^Hy\ 

.f-I^^.Sl (-2) 



9x \ ddxdy dx"^ 
Different physical situations arise for different boundary conditions. We shall con- 
sider first the half plane problem bounded by the plane y = 



5 Reflection of elastic waves from a plane boundary 

Consider the half space y > 0, — oo < x < oo. Several types of boundary conditions 
can be prescribed on the plane boundary : (i) dynamic: the stress components only 
(the traction condition); (ii) kinematic: the displacement components only, or (iii). a 
combination of stress components and displacement components. Most difficult are 
(iv) the mixed conditions in which stresses are given over part of the boundary and 
displacements over the other. 

We consider the simplest case where the plane y = is completely free of external 
stresses, 

'^yy ~ — 0, (5.23) 

and 

V = (5.24) 
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It is clear that (5.23) affects the P and SV waves only, while (5.24) affects the SH 
wave only. Therefore we have two uncoupled problems each of which can be treated 
separately. 

5.1 P and SV waves 

Consider the case where only P and SV waves are present, then — Hy — 0. Let all 
waves have wavenumber vectors inclined in the positve x direction: 

= f{y)e'^''~'^\ = hMe"^^-'^' (5.25) 
It follows from (4.15) and (4.16) that 

|/^„V^0, ^+/»%^0, (5.26) 

with 



o^=\ hr-e-x kl-e, P=\ -^-e^\ k'T-C' (5.27) 



'-L V T 



We first take the square roots to be real; the general solution to (5.26) are sinusoids, 
hence. 

On the right-hand sides the first terms are the incident waves and the second are the 
reflected waves. If the incident amplitudes Ap, As and are given, what are the properties 
of the reflected waves Sp, S5? The wave number components can be written in the polar 
form: 

(^,a) = A;L(sin^L,cos^L), ((, /9) = A;T(sin^T, cos^t) (5.29) 

where (/cl, kx) are the wavenumbers, the {9l, Ot) the directions of the P wave and SV 
wave, respectively. In terms of these we rewrite (5.28) 

_ j^^^/tkT(sm.9TX—cos6Ly—'^t) _j_ ^^ikT{sm6TX+cos6Ty—0Jt) g-j^^ 

In order to satisfy (5.23) {Tyy — T^y = 0) on y = for all x, we must insist: 

kL sin 9l = kr sin 9t, (^ = C) (5-32) 



3.5. Reflection of elastic waves from a plane boundary 
This is in the form of Snell's law: 

sin 01 sin 9t 
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cl 



Ct 



implying 



sin6'L Cl / A + 2/x kj 



K 



(5.33) 



(5.34) 



sin^T Ct \l IJ' ki 
When (5.23) are applied on y — Q the exponential factors cancel, and we get two 
algebraic conditions for the two unknown amplitudes of the reflected waves {Bp, Bs) : 



A;i(2sin2 Ol - K^){Ap + Bp) - 4sin2^T(^5 - S5) = 
kl sin 2^i(Ap - Bp) - k^ cos 9t{As + Bs) = 0. 

Using (5.34), we get 

2 sin^ Ol-k^ = K^{2 sin^ Or - I) ^ -k^ cos 2^t 
The two equations can be solved and the solution expressed in matrix form: 



(5.35) 
(5.36) 





V 


Spp 


Ssp 






Sps 


Sss 




(5.37) 



S = 



where 

Spp Ssp 
Sps Sss 

denotes the scattering matrrix. Thus Sps represents the reflected 5"- wave due 
P wave of unit amphtude, etc. It is straightforward to verify that 

sin 201 sin 2^t - k,^ cos^ 2^t 



Spp — 



•>pp 

Ssp 
Sps 
Sss 



sin 261 sin 2^t + cos^ 2^t 

-2/«2 sin 29t cos 29t 
sin 201 sin 29t + cos^ 26't 

2 sin 26Il cos 26It 
sin 29 L sin 26*^ + k,^ cos^ 26't 
sin 2^L sin 29t - cos^ 29t 
sin 2^L sin 2^t + cos^ 2^t 



(5.38) 
to incident 

(5.39) 

(5.40) 
(5.41) 
(5.42) 
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In view of (5.33) and 



The scattering matrix is a function of Poisson's ratio and the angle of incidence. 

(i) P- wave Incidence : In this case is the incidence angle. Consider the special 
case when the only incident wave is a P wave. Then Ap and As — and only 
Spp and Ssp are relevant. . Note first that > in general . For normal incidence, 
Ol = 0, hence Or = 0. We find 

Spp = -1, Sps = (5.44) 

there is no SV wave. The refelcted wave is a P wave. On the other hand if 

sin 29l sin 29t - cos^ 29t = (5.45) 

then Spp = 0, hence Bp = but Bs 7^ 0; only SV wave is reflected. This is the case 
of mode conversion, whereby an incident P waves changes to a SV wave after reflection. 
The amplitude of the reflected SV wave is 

j; = Sps= (5.46) 



(ii) SV wave Incidence : Let Ap = but As ^ 0. In this case 6t is the incidence 
angle. Then only Ssp and Sss ai'e relevant. For normal incidence, 6l = 6t = 0, 
Sss — and Ssp = 0; no P wave is reflected. Mode conversion {Bp ^ 0, Bs = 0) also 
happens when (5.45) is satisfled. Since 9l > 9t, there is a critical incidence angle 9t 
beyond which the P wave cannot be reflected back into the solid and propagates only 
along the x axis. At the critical angle 

sin^L = 1, (^L = 7r/2), or sin^r = l/K (5.47) 

by Snell's law. Thus for u = 1/3, k = 2 and the critical incidence angle is 9t = 30°. 

Beyond the critical angle of incidence, the P waves decay exponentially away from 
the free surface. The amplitude of the SV wave is linear in y which is unphysical, 
suggesting the limitation of unbounded space assumption. 
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Figure 4: Amplitude ratios for incident P waves for various Possion's ratios. From Graff: 
Waves in Elastic Solids. Symbols should be converted according to : Ai ^ Ap, A2 ^ 
Bp, Bi As, B2 Bs. 
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Figure 5: Reflected wave amplitude ratios for incident SV waves for various Possion's 
ratios. Prom Graff: Waves in Elastic Solids. Symbols should be converted according to 

: Ai — > Ap, A2 — > Bp, Bi — > As, B2 Bs- 
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5.2 SH wave 

Because of (5.2) 



dH^ dHy ^ 
+ = 



dx dy 

we can introduce a stream function ip so that 



where 



Clearly the out-of-plane dispacement is 



and 



^ dy dx dx"^ dy'^ 

r,. = /.^V^ = ^-^ (5.51) 



The zero-stress boundary condition implies 

Thus the problem for ip is analogous to one for sound waves reflected by a solid plane. 
Again for monochromatic incident waves, the solution is easily shown to be 

i/j = {Ae-'f^y - Ae'f^y) e^°=^-^'^* (5.53) 

where 

We remark that when the boundary is any cylindrical surface with axis parallel to 
the z axis, the the stress-free condition reads 

Tzn — 0, on B. (5.55) 

where n is the unit outward normal to B. Since in the pure SH wave problem 

^ _ ^ _ Ay^^ - — — ^ 
^"^ ^ dn dn dn dt'^ 

Condition (5.55) implies 

^ = 0, on B. (5.56) 

on 

Thus the analogy to acoustic scattering by a hard object is true irrespective of the 
geomntry of the scatterer. 
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6 Rayleigh surface waves 

Refs. GrcLfF, Achenbach, Fung 

In a homogeneous elastic half plane, in addition to P, SV and SH waves, another 
wave which is trapped along the surface of a half plane can also be present. Because 
most of the action is near the surface, this surface wave is of special importance to 
seismic effects on the ground surface. 

Let us start from the governing equations again 

dx^ dy^ cldt^' ^ ' ' 

We now seek waves propagating along the x direction 

= (/(y)e^€^-''^*) , H, = ^{h{y)e'^''-'''*) (6.3) 
Then f{y),h{y) must satisfy 

+ (-V4-f)/ = o, + (.v4-«=)fc = o, (6.4) 

To have surface waves we insist that 



a^^e- c^v4, p-^e- u^yc^T (6.5) 

be real and postive. Keeping only the solutions which are bounded for y oo, we get 

(/) = ^e-"f e^(^^-'^*\ //^ = Se-^V(^^-'^*l (6.6) 

The expressions for the displacements and stresses can be found straightforwardly. 

= (^i^Ae-^'y -pBe-^y^e'^^''-^'\ (6.7) 

Uy = - (aAe-^y + iCBe-^y^ e'^^''-'''\ (6.8) 

rxx = fi{(^p^ -^^ -2a^^Ae-''y -2iP^Be-^y^e'^^''-^^\ (6.9) 

^yy = + C^^Ae-^y + 20CBe-^y^e'^^''-''*\ (6.10) 

T^y = /X |-2meyle-"^ + + P^) Be-f^y^ e^(«^-'-*) (6.11) 
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On the free surface the traction-free conditions Tyy — Txy — require that 

(p' + e)A + 2tP^B = 0, (6.12) 
-2iaCA + (p^ + B = 0. (6.13) 

For nontrivial solutions of A, B the coefficient determinant must vanish, 

(^' + f)'-4a:^f = 0, (6.14) 



or 

,2 



/-2 



Cfj-\ 




(6.15) 



which is the dispersion relation between frequency cu and wavenumber ^. Prom either 
(6.12) or (6.13) we get the amplitude ratio: 



B ^2 2ia^ 

In terms of the wave velocity c — cu/^, (6.15) becomes 



(6.16) 



or, upon squaring both sides, finally 




(6.18) 



where 



cl \ + 2a 2-2u 

k — 



ct y V 1 — 2i/ 

The first solution c = o; = is at best a static problem. In fact a = (3 = ^ and A = —iB, 
so that Ux — Uy = which is of no interest. 

We need only consider the cubic equation for c^. Note that the roots of the cubic 
equation depend only on Poisson's ratio, through — 2(1 — — 2v). There can be 
three real roots for c or a;, or one real root and two complex-conjugate roots. We rule 
out the latter because the complex roots imply either temporal damping or instability; 
neither of which is a propagating wave. When all three roots are real we must pick the 
one so that both a and /? are real. We shall denote the speed of Rayleigh wave by cr. 
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Figure by MIT OCW. 

Figure 6: The velocity of Rayleigh surface waves cr. From Fung Foundations of Solid 
Mechanics. 



For c = 0, the factor in curley brackets is 



{.} = -16 1- 



< 



For c = ct the same factor is equal to unity and hence positive. There must be a solution 
for c such that < c < ct- Furthermore, we cannot have roots in the range c/ct > 1- 
If so, 



< 



which is not a surface wave. Thus the surface wave, if it exists, is slower than the shear 
wave. 

Numerical studies for the entire range of Poisson's ratio (0 < z/ < 0.5) have shown 
that there are one real and two complex conjugate roots if u > 0.263 . . . and three real 
roots if u < 0.263 . . .. But there is only one real root that gives the surface wave velocity 
cr. a graph of cr for all values of Poisson's ratio, due to Knopoff , is shown in Fig. 6. 
A curve-fitted expression for the Rayleigh wave velocity is 



Cij/CT = (0-87 + 1 ■12z/)/(l + z/) 



(6.19) 



For rocks, A = /i and = \, the roots are 

{c/crf = 4, 2 + 2/V3, 2 - 2/V3. 



(6.20) 
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The only acceptable root for Rayleigh wave speed cr is 

{cR/crf = (2 - 2/V3)^ = • 9194 



(6.21) 



or 



Or = 0.9588CT. (6.22) 
The particle displacement of a particle on the free surface is, from (6.7) and (6.8) 



2e 



Note that 



A 



p +e 

2^ 



A 



> 



b^A 



—a + 



p +e 

2e 



A 



2e 



> 



hence 



and 



Ux — asin(a;i — ^x), Uy — hcos{ujt — ^x) 



+ = 1 

a2 ^ 62 



The particle trajectory is an ellipse. In complex form we have 



^ + = exp {i {ujt -^x- n/2)} 
a b 



(6.23) 



(6.24) 



(6.25) 



(6.26) 



Hence as t increases, a particle at {x, 0) traces the ellipse in the counter-clockwise direc- 
tion. See figure (6). 



3. 7. Scattering of SH waves by a circular cavity 



24 



Instantaneous 

Instantaneous wave surface Direction of wave propagation 
particle velocity \ >- 




Particle path 



Figure by MIT OCW. 

Figure 7: Displacement of particles on the ground surface in Rayleigh surface wave From 
Fung Foundations of Solid Mechanics. 

7 Scattering of monochromatic SH waves by a cav- 
ity 

7.1 The boundary- value problem 

We consider the scattering of two-dimensional SH waves of single frequency. The time- 
dependent potential can be wirtten as 



where the potential is governed by the Helmholtz equation 



CO 



(7.1) 



(7.2) 



To be specific consider the scatterer to be a finite cavity of some general geometry. On 
the stress-free boundary B the shear stress vanishes, 



-3? 



''^ ^ on 



-iujt 







hence 



dn 



0, on B 



Let the incident wave be a plane wave 

0j = Ae'*^-" 



(7.3) 

(7.4) 
(7.5) 
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Let the angle of incidence with respect to the positive x axis be 9o- In polar coordinates 
the incident wave vector is 

k = A;(cos^o, sin^o), x = r(cos^,sin^) (7.6) 

and the potential is 

(f)i^A exp [iA;r(cos Oo cos 9 + sin Oo sin 6)] = Ae^'=''"°^(^-^°^ (7.7) 
Let the total wave be the sum of the incident and scattered waves 

(p = (pi + (ps (7.8) 
For the scattered wave the boundary condition on the cavity surface is 

dn dn ^ ' 
In addition, the scattered wave must satisfy the radiation condition at infinity, i.e., it 
must propagate outward at infinty. 

Note that this boundary value problem is identical to the one for sound plane sound 
wave in air (or water) scattered by a perfectly rigid cylinder. 

We treat below the special case of a circular cavity of readius a This is one of the 
few geomentries that can be solved analytically. 



7.2 The circular cavity 

In polar coordinates the governing equation reads 

r dr \ dr J dO"^ ' 
Since 0/ satisfies the preceding equation, so does 05 . 

First, It is shown in Appendix A that the plane wave can be expanded in Fourier- 
Bessel series : 

oo 

^ikrco.{e-eo) ^ J2 Sni^Jnikr) COS n(e - Oo) (7.11) 

n=0 

where e„ is the Jacobi symbol: 

eo = 0, e„ = 2, n= 1,2,3,... (7.12) 
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Each term in the series (7.11) is called a partial wave. 
By the method of separation of variables, 

Mr,o) = R{r)e{e) 

we find 

r^R" + rR' + {k\^ -n'^)R = 0, and 9" + n^e = 

where n = 0, 1, 2, . . . are eigenvalues in order that 6 is periodic in 9 with period 27r. 
For each eigenvalue n the possible solutions are 

Qn = {sin n9 , COS n6) , 

Rr.^{H^'\kr),Hi'\kr)), 

where Hn\kr) , Hn\kr) are Hankel functions of the first and second kind, related to 
the Bessel and Weber functions by 

Hi'\kr) = Jn{kr) + iY^ikr), H^^\kT) = J„(A;r) - lYnikr) (7.13) 

The most general solution to the Helmholtz equation is 

oo 

05 = A J](A„sinne + S„cosne) [C„//(i)(A;r) + L>„//f (A;r)] , (7.14) 

n=0 

For large radius the asymptotic form of the Hankel functions are 

" V Tr/cr " V Tr/cr 

In conjunction with the time factor exp(— icjt), Hn^ gives an outgoing wave while h!^'^ 
gives an incoming wave. To satisfy the radiation condition, we must discard all terms 
involving Hn'' ■ From here on we shall abbreviate Hn'^ simply by H^- The scattered 
wave is now 

oo 

05 = ^ ^ {An sin nO + cos nO) Hn{kr) (7.16) 

n=0 

The expansion coefficients {An, Bn) must be chosen to satisfy the boundary condition 

on the cavity surface^ Once they are determined, the wave is found everywhere. In 

"'^In one of the numerical solution techniques, one divides the physical region by a circle enclosing the 
cavity. Between the cavity and the circle, finite elements are used. Outside the circle, (7.16) is used. 
By constructing a suitable variational principle, finite element computation yields the nodal coefficients 
as well as the expansion coefficients. See (Chen & Mei , 1974). 
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particular in the far field, we can use the asymptotic formula to get 



05 ~ >1 ^ {An sin n9 + Bn cos n9) e' 



■inn/ 2 



n=Q 



^ikr—iir/A: 



Tikr 



(7.17) 



Let us define the dimensionless directivity factor 

oo 

A{e) = {An sin ne + S„ cos nO) e'^'^'^l'^ 



(7.18) 



n=0 



which indicates the angular variation of the far-field amplitude, then 



05 ~ AA{Q) 



J.ki — i7r/4 



7r/cr 



(7.19) 



This expression exhibits clearly the asymptotic behaviour of 05 as an outgoing wave. 
By differentiation, we readily see that 



lim \/r 



dr 







(7.20) 



which is one way of stating the radiation condition for two dimensional SH waves. 
Let us complete the solution. 

Without loss of generahty we can take Qo — ^- On the surface of the cylindrical 
cavity r = a, we impose 



90/ <905 „ 

or or 



It follows that A„ — and 



e''i''AJ'^{ka) + BnkH'^{ka) = 0, n = 0, 1, 2, 3, . . . n 



where primes denote differentiation with respect to the argument. Hence 



B - -Af '^»(^") 



H'n{ka) 



The sum of incident and scattered waves is 

oo 

= ^ ^ e„ 

and 



n=0 



Jn{kr) - ^^Mkr) 



H'n{ka) 



cos n9 



n=0 



cos n9 



(7.21) 



(7.22) 
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Figure 8: Angular distribution of scattered energy in the far field in cylindrical scattering 

The limit of long waves can be approximatedly analyzed by using the expansions for 
Bessel functions for small argument 

2 n 

Jo{x) ~ 1 - J + 0{x% Jn{x) - ;n = 1, 2, 3... 
Y^{x)^-\ogx, Y^{x) ^^^,71 = 1,2,3... (7.23) 

Then the scattered wave has the potential 

= ^{kaf (^-^-Ho{kr) - Hi{kr) cosO^ + 0{kaf (7.24) 

The term HQ{kr) coresponds to a oscillating source which sends istropic waves in all 
directions. The second term is a dipole sending scattered waves mostly in forward and 
backward directions. For large kr, the angular variation is a lot more complex. The far 
field pattern for various ka is shown in fig 4. 

On the cavity surface surface, the displacement is proportional to V'^tp{a, 9) or 
V^0(a, 9). The angular variation is plotted for several ka in figure 5. 

For numerical simulations, see the website 
http: / /ocw.mit.edu/OcwWeb/Civil-and-Environmental-Engineering/l-138JFall-2004/Simulations/ 
Remark on energy flux: At any radius r the total rate of energy outflux by the 
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Figure 9: Polar distribution of 0(a, ^) on a circular cylinder. 
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scattered wave is 



i-2n 

7. 



dt 

/looker 



p2tt 




Jo 


or 



fjwk'^r 

~ z <J 



U [iuk'^(t)e-^^] 
dO (t)* ■ 

or 



(7.25) 



where overline indicates time averaging over a wave period 2t:/uj. The integral can be 
evaluated by using the asymptotic expression. 

Remark: In the analogous case of plane acoustics where the sound pressure and 
radial fluid velocity are respectively, 



(90 , dcj) 
P = ~Po-Q^^ and 



(7.26) 



the energy scattering rate is 

f27r 



^0 



dOpUr 



3? / d9[-i 



c 



dr 



^ de((f)^ 



c 



dr 



(7.27) 



8 The optical theorem 

For the same scattcrcr and the same frequency uj, different angles of incidence 6j define 
different scattering problems 4>j. In particular at infinty, we have 



Let us apply Green's formula to 4>i and 02 over a closed area bounded by a closed 
contour C, 

//^(^.vV,-*-v^*-.)<(^ = /^(0.^-.^,^)d. + {^^jt 

where n refers to the unit normal vector pointing out of 5". The surface integral vanishes 
on account of the Helmholtz equation, while the line integral along the cavity surface 
vanishes by virture of the boundary condition, hence 

(901 (902" 



'2 



dn dn 



ds^O (8.2) 
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By similar reasioning, we get 



ds 



^ dn ^ dn 
where 0* denotes the complex conjugate of 0i. 
Let us choose 0i = 02 = in (8.3), and get 



ds = 



(8.3) 



ds 



c 



ds = 2» 



ds ( 



dn 







Physically, across any circle the net rate of energy flux vanishes, i.e., the scattered power 
must be balanced by the incident power. 
Making use of (8.1) we get 



c2tt 

= 9/ rde 



I 

Jo 



^ikrcos(9-9o) ^ / ^^(^)e^*^^-^'^/^ 



nkr 



-ikcos(9 - 0\e-^''''''"^^^-'^°^ - ik\ —Al{e)e-'^'+'^'^ 

nkr 



= 9 / rde l-ikcos(e - 9o) + -^{-ik)\Ao\ 
Jo { TT/cr 



The first term in the integrand gives no contribution to the integral above because of 
periodicity. Since Q{if) — Q{if*), we get 







— / \Ao{9)fd9 

rde !^Aoi-ik)^l^[l + cos{e - ^^)]eW4gifcr(i-cos(e-e<,)) 



271 



Ao(e)\^de 



For large kr the remaining integral can be found approximately by the method of sta- 
tionary phase (see Appendix B), with the result 



/■27r 

/ deAo{0)[l + cos{e - 0^)]e"'r{i-cos{e-e„)) ^ ^^^^^ 
Jo 



27r 




iTv/4 



(8.5) 
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We get finally 



(8.6) 



Thus the total scattered energy in all directions is related to the amplitude of the 
scattered wave in the forward direction. In atomic physics, where this theorem was 
originated (by Niels Bohr) , measurement of the scattering amplitude in all directions is 
not easy. This theorem suggests an econmical alternative. 

Homework For the same scatterer, consider two scattering problems 0i and 02- 
Show that 



For general elastic waves, see Mei (1978) for similar and other identities in elastody- 
namics with rigid inclusions. J . Acoust. Soc. Am. 64(5), 1514-1522. 

9 Diffraction of SH wave by a long crack - the parabolic 
approximation 

References 

Morse & Ingard, Theoretical Acoustics Series expansions. 

Born & Wolf, Principle of Optics Fourier Transform and the method of steepest de- 
scent. 

B. Noble. The Wiener-Hopf Technique. 

If the obstacle is large, there is always a shadow behind where the incident wave 
cannot penetrate deeply. The phenomenon of scattering by large obstacles is usually 
referred to as diffraction. 

Diffraction of plane incident SH waves by a long crack is identical to that of a hard 
screen in acoustics. The exact solution was due to A. Sommerfeld. We shall apply 
the boundary layer idea and give the approximate solution valid far away from the tip 
kr ^ 1 hy the parabolic approximation, due to V. Fock. 

Refering to figure () let us make a crude division of the entire field into the illuminated 
zone I , dominated by the incident wave alone, the reflection zone II dominated the sum 
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Figure 10: Wave zones near a long crack 

of the incident and the reflected wave, and the shadow zone III where there is no wave. 
The boundaries of these zones are the rays touching the crack tip. According to this 
crude picture the solution is 

Af, exp{ik cos 9x + ik sin 9y) , I 
Ao[cxp{ik cos 9x + ik sin 6y) + exp{ik cos 9x — ik sin 9y)], II (9.1) 
0, /// 

Clearly (9.1) is inadquatc because the potential cannot be discontinuous across the 
boundaries. A remedy to provide smooth transitions is needed. 

Consider the shadow boundary Ox'. Let us introduce a new cartesian coordinate 
system so that x' axis is along, while the y' axis is normal to, the shadow boundary. 
The relations between (x, y) and (x', y') are 

x' — xcos9 + ysin9, y' — y cos9 — xsin9 (9.2) 

Thus the incident wave is simply 

07 = A„e^*^^' (9.3) 
Following the chain rule of differentiation, 

90 (90 dx' d4> dy' ^90 ■ 

dx dx' dx dy' dx dx' dy' 
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d(j) d(j) dx' d(t) dy' ■ q^^ 
dy dx' dy dy' dy dx' dy' 

we can show straightforwardly that 

d^(f) d^(j) _ d^(f) d'^ct) 
dx"^ dy"^ dx'"^ dy'"^ 

so that the Helmholtz equation is unchanged in form in the x', y' system. 

Wc try to fit a boundary layer along the x' axis and expect the potential to be almost 
like a plane wave 

4>{x: ,y')^A{x^,y'y^-' (9.4) 

, but the amplitude is slowly modulated in both x' and y' directions. Substituting (9.4 
into the Helmholtz equation, we get 

Expecting that the characteristic scale of A along x' is much longer than a wavelength, 
kLj. ^ 1, we have 

dA ^ ^ ^ dA d^A 

— < kA, hence 2z/c— > — 
dx dx dx^ 

We get as the first approximation the Schrodinger equation^ 

, dA d^A , , 

In this transition zone where the remaining terms are of comparable importance, hence 
the length scales must be related by 

k 1 

x' y'' 

Thus the transition zone is the interior of a parabola. 

Equation (9.7) is of the parabolic type. The boundary conditions are 



^2, implying ky' Vkx' 



A(x,oo)=0 (9.8) 



^In one-dimensional quantum mechanics the wave function in a potential-free field is governed by 
the Schrodinger equation 

., dtp 1 „ 
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A{x, -oo) = Ao (9.9) 

The initial condition is 

AM) J (9.10) 

[ ^0, y' < 

he initial-boundary value for A has no intrinsic length scales except x', y' themselves. 
Therefore the condition kL^ » 1 means kx' » 1 i.e., far away from the tip. This 
problem is somwhat analogous to the problem of one-dimensional heat diffusion across 
a boundary. A convenient way of solution is the method of similarity. 
Assume the solution 

A = AJi^) (9.11) 

where 

yirkx 

is the similarity variable. We find upon subsitution that / satisfies the ordinary differ- 
ential equation 

f"-i7rjf'^0 (9.13) 
subject to the boundary conditions that 

/^O, 7^-oo; /^l, 7^oo. (9.14) 

Rewriting (9.13) as 



f" 

we get 



log/' = i7i'j/2 + constant. 

One more integration gives 



f = C I exp ) 



Since 



au 



2 J V2 



we get 



g-i7r/4 
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and 
/ 



A 



exp 



du 



Defining the cosine and sine Fresnel integrals by 



C{l) 
we can then write 



cos 



TTV 



dv, 5(7) 



sm 



exp 



71V 



du } (9.15) 



dv 



V2 



+ i 



+ 5(7) 



(9.16) 



(9.17) 



In the complex plane the plot of C{'y) +iS{'y) vs. 7 is the famous Cornu's spiral, shown 
in figure (9). 

The wave intensity is given by 



+ C(7) 



+ 



+ 5(7) 



(9.18) 



Since C,S — > —1/2 as 7 — > —00, the wave intensity diminshes to zero gradually into 
the shadow. However, C, 5" 1/2 as 7 ^ 00 in an oscillatory manner. The wave 
intensity oscillates while approaching to unity asymptotically. In optics this shows up 
as alternately light and dark diffraction bands. 

In more complex propagation problems, the parabolic approximation can simplify the 
numerical task in that an elliptic boundary value problem involving an infinite domain 
is reduced to an initial boundary value problem. One can use Crank-Nicholson scheme 
to march in "time", i.e., x'. 

Homework: Find by the parabolic approximation the transition solution along the 
edge of the reflection zone. 
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Figure 11: Cornu's spiral, a plot of the Fresnel integrals in the comiplex plane of C(7) + 
iS{j). Abscissa: C(7). Ordinate : iS^j). 
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Figure 12: Diffraction of a normally incident plane sound wave on a rigid half screen. 
Ordinate: /(7) = A/Aq. Abscissa: —7 = ky' / \/ nkx' . 
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10 Exact theory of Wedge Diffraction 

Refs. J. J. Stoker 1957 , Water waves, pp 120-125 
Born & Wolf, 1950 Principles of Optics 
Noble, 1963 The Wiener-Hopf Technique 

The diffraction of plane incident waves by a semi-infinite barrier is a celebrated 
problem in classical physics and was first solved by A. Sommerfeld. Several analytical 
treatments are available, including the mathematically very elegant technique of Wiener 
& Hopf is available . In the last section we gave an approximate theory by parabolic 
approximation. Here we present an exact theory for the more general case of a wedge, 
by the more elementary method of series expansion. Though the result is in an infinite 
series and not in closed form, quantitative information can be calculated quite readily 
by a computer. 

The basic ideas were described for water waves by Stoker (1957) for a wedge. Ex- 
tensive numerical results were reported by Dr. H.S. Chen (Army Corps of Engineers, 
Tech Rept: CERC 878-16, 1987). Based on Stoker's analysis more extensive numerical 
computations have been carried out by Dr. CD. Li. To facilitate the understanding 
of the physics, these results are presented in animated form in our subject website (see 
Simulations) : 

http: / /ocw.mit.edu/OcwWeb/Civil-and-Environmental-Engineering/l-138JFall-2004/ 

We shall describe the problem for water waves; the analysis and results can be 
adapted for cracks, sound,... etc. 

Refering to Figure 10, we consider a vertical wedge of arbitrary apex angle in a 
sea of constant depth h. Let the tip of the wedge be the z axis orgin and the x axis 
coincide with one wall. The still water surface is the x, y plane. In the cylindrical polar 
coordinate system (r, 9, z) with {x = r cos 9,y = r sin 9), z, the walls are given by ^ = 
and i/TT with 1 < u < 2. A train of monochromatic waves is incident from infinity at the 
angle a with respect the x axis. 

In the water region defined hy < 9 < un and > 2; > —h, the velocity potential 
$(r, 9, z, t) must satisfy the Laplace equation. 



1 9$ 



1 



(10.1) 



dr'^ r dr 



+ 



,2 g02 g^2 



= 



r 
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X 



Figure 13: Coordinate system and subregions. 



and subject to the linearized free surface boundary conditions 



9$ 

~dt 



dC _ <9$ 
dt dz 



whicdh can be combined to give 



dz 



0, z = 0. 



(10.2) 
(10.3) 

(10.4) 



Along the impermeable bottom and walls, the no flux boundary conditions are 

9$ 



dz 







9$ 

The incident wave train is given by 



on z = —h 



at 9 = and z/tt 



-igAo cosh k{z + h) 



r, y e 



-ikr cos{9—a)—iu)t 



ijj cosh kh - - 
where k is the real wavenumber satisfying the dispersion relation 



(10.5) 
(10.6) 

(10.7) 



tu^ = gk tanh kh, 



(10.8) 
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and 77 + a is the angle of incidence measured from the x axis. Aq is the incident wave 
amplitude. 

Because of the vertical side-walls, the three dimensional problem can be reduced to 
a two dimensional one by letting 

where 0(r, 6, t) is the horizontal pattern of the velocity potential normlized for an inci- 
dent wave of unit amplitude. 

Substituting Equation (10.9) into the Laplace equation and using both the kinematic 
and dynamic boundary conditions on the free surface, the Laplace equation is then 
reduced to the Helmholtz equation 

g^ + ^^ + ^a| + «^^* = 0. -S»SO. (10.10) 
with the following boundary conditions on the rigid walls of the wedge: 

^ = at 6 = and 00 (10.11) 

09 

The free surface displacement C from the mean water level z=0 can be represented by 

C(r,^,t) = -^^=Aov{r,e)e-'^' (10.12) 
9 ot 

Note that 77 is dimensionless. 

Refering to Figure 10, the entire water region can be divided into three zones accord- 
ing to the crude picture of geometrical optics. 1: the zone of incident and reflected plane 
waves, II : the zone of incident plane wave and III ; the shadow with no plane wave. 
This crude picture is discontinuous at the border lines separating the zones. Siknce the 
physical solution must be smooth every where we must find the transitions. Let us use 
the ideas of boundary layers. 

Let the total potential be expressed in a compact form by 

= 0^(r, 9) + 0,(r, 9), for all < 9 < i^n (10.13) 

where 0o consists of only the plane waves, 

(f)i + (f)r n — a > 9 > 0, in I; 
Mr,0)={ (Pi Tr + a>9>7r-a, mil; (10.14) 

9o > 9 > IT + a, in III. 
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Here 0j denotes the incident wave 



(f>i = e 



-ikr cos{&— a) 



where a denotes the angle of incidence, and 4>r the reflected wave 



(f)r = e 



—ikr cos{0+a) 



(10.15) 



(10.16) 



The correction is the diffracted wave 0^ which must satisfy the radiation condition and 
behaves as an outgoing wave at infinity, i.e.. 



or 



hm A/r(^ - ik(l)s) = 
r^o Or 

r — > oo 



kr 



(10.17) 



(10.18) 



10.1 Solution by Fourier series (or, finite Fourier Transform) 

Let us solve the scattered wave formally by Fourier series 



1 - 2 °° - 

0(r,^) = — 0o(r) + — V4 



nd 



r cos 



;i0.19) 



n=l 



then, the Fourier coefficients are: 

(pnikr) 



nd 

(f){kr, 9) cos — d9 



V 



;io.2o) 



From (10.10), each Fourier coefficient satisfies 



h r h 



The general solution finite at the origin is 



(t)n{kr) = anJn/u{kr) 

where the coefficient's On, n = 0, 1, 2, 3, ... are to be determined. 
The Fourier coefficient of (10.13) reads 

4>n{kr) — anJn/i,{kr) — / ^^cos — d9 + / ^qCOS — d9 

Jo ^ Jo ^ 



(10.21) 



;i0.22) 



(10.23) 
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or 



(ps,n = anJn/uikr) - (f)o,n (10.24) 

Applying the operator liiUr^oo V^{d /dr — ik) to both sides of (10.23), and using the 
Sommerfeld radiation condition (10.17), we have 



hm \/r ( 

1 — ►oo V 



f d_ 
dr 



— ik 



, COS — ad 







;i0.25) 



We now perform some asymptotic analysis for large kr to evaluate 
The first part can be treated explicitly for large kr, since 



Jn/u{kr) 



nkr 



, , nn 7r\ 

cos kr 

2u 4/ 



It follows that 



lim \/r 



d_ 
dr 



ik ) Jn/u{kr) 




2^e-<fer-|^+f) 

TT 



;i0.26) 
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For the second part, we substitute 0o from (10.15) and (10.16) to rewrite the integral as 

1 2 



■ COS — dQ 



V 



^0 



e-^'='-^°«('^-")cos— + 



g-ifcr-cos(e+a) cos— 
1/ 



+ / e-^'^'-^^^^^-^^cos— 

V 



;i0.28) 



Each of the integrals above can be evaluated for large kr by the method of stationary 
phase (again). Details are given in the next subsection, only the results are cited below. 
The first integral is appoximately 



2tt' 
kr 



from which 



lim i/r [ — ik 



kr 
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;i0.29) 



g-ifer-cos(e-a) cos—de 



lim y/r 



d_ 
dr 



— ik] I cos 
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2v^cos(^)e-'^'- 



(10.30) 
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where we have used i — e^'^^^. By similar analysis the second integral is found to be 



12(0) - cos 



It follows that 
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(10.31) 



;i0.32) 



Finally the third integral is approximately 



hence 
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(10.33) 



(10.34) 



In summary, only the first integral associated with the incident wave furnishes a 
nonvanishing contribution to the expansion coefficients, i.e., 

lim f I- - ik] r 0„cos -de ~ 2v^cos !^e-<'='^+^) (10.35) 
r--»oo \or J Jq V V 

With this result we get by substituting (10.27) and (10.35) into (10.25), the coefficients 
a„ are found 

a„ = 27r cos — (10.36) 
By the inverse transform, (10.19), we get the exact solution, 

2 



0(r,^) 



V 



Jo[kr) + 2 2^ e 2.. j^i^[kr) cos — cos 

n=l ^ 



(10.37) 



Numerical computations by MATLAB is straightforward. Aside from its own practical 
interest, this solution is useful for checking strictly numerical methods for diffraction 
problem for other geometries. 



3.10. Exact theory of Wedge Diffraction 



44 



10.2 Asymptoic approximation of Integrals 

For the first integral Ji, we take the phase to be fi{0) = kcos{9 — a). The points of 
stationary phase must be found from 

f[(e) = -k sm(e - a) = 0, (10.1) 

hence 9 = a,a ±7t. Only the first at 6*1 = a lies in the range of integration (0, tt — a) 
and is the stationary point. Since 

f'^(e^) = -k cos(ei -a)^-k<0 (10.2) 

the integral is approximately 



/^(^)~cos(^)e-^'='-^°^(^^-")+^ 
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kr 



cos( — )e 



-ikr+- 
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(10.3) 



For the second integral /2, we take the phase to be f2{d) — kcos{9 + a). The 
stationary phase point must be the root of 

/^(^) = -A;sin(^ + a) = (10.4) 

or ^ = —a, —a ± i. The stationary point is at ^2 = tt — a which is the the upper limit 
of integration. Since 

f!^(92) = -k cos(e2 + a)^k>0 (10.5) 

I2 is approximately 



/2(^)~^cos(^)e-^^^^°^(^^+")-'^ 
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(10.6) 



Lastly for the third integral /3, the phase is fz{9) — kcos{9 — a). The point of 
stationary phase is found from 

f'^{9) = -k sm(9 - a) = (10.7) 

or ^ = TT, ±7r + a. Only the point ^3 = tt + a is acceptable and coincides with the upper 
limit of integration. Since 

/^(^s) = -A;cos(^3 -a)^k>0, (10.8) 

73 is approximately 



13(9) ^ - cos( 



kr 
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10.3 Two limiting cases 

(1) A thin barrier. Let the wedge angle be by setting u = 2. Equation (10.37) then 
becomes 

oo „ 

0(r, e) = Jo{kr) + 2 J2 e""^ Jn/2{kr) cos ^ cos y (10.10) 

n=l 

(see Stoker (1957)). 

(2) An infinite wall extending from x = —oo to oo. Water occupying only the half 
plane of y > and the wedge angle is 180 degrees. The diffracted wave is absent from 
the solution, and the total wave is only the sum of the incident and reflected waves : 

(/)(r, e) = e-^'=''^°'('^-") + e-^^'-'=°^(^+") (10.11) 

By employing the partial- wave expansion theorm, (Abramowitz and Stegun 1964), the 
preceding equaion becomes 



;io.i2) 



Jo{kr) + 2 ^^(— i)" J„(A;r) cos no; cos 

n=l 

which agree with (10.37) for u = 1. 

For the animated version of the wave patterns for different wedge angles snd angle 
of incidence, please visit the website. 
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A Partial wave expansion 



A useful result in wave theory is the expansion of the plane wave in a Fourier series 
of the polar angle 9. In polar coordinates the spatial factor of a plane wave of unit 
amplitude is 



Consider the following product of exponential functions 



^zt/2^-z/2t 
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n\ l!(n + l)! 2!(n + 2)! 
The coefficient of is nothing but Jniz), hence 
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r\(n + r)\ 



exp 
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t 



Now we set 



t = ie z = kr. 



The plane wave then becomes 

oo 
N=-oo 

Using the fact that J_„ = (— l)"Jn, we finally get 

oo 



n=0 



(A.l) 



where e„ is the Jacobi symbol. The above result may be viewed as the Fourier expan- 
sion of the plane wave with Bessel functions being the expansion coefficients. In wave 
propagation theories, each term in the series represents a distinct angular variation and 
is called a partial wave. 

Using the orthogonahty of cos nO, we may evaluate the Fourier coefficient 



Jn{kr) = — 



Akr cos 9 



COS nOdO, 



(A.2) 



which is one of a host of integral representations of Bessel functions. 
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B Approximate evaluation of an integral 

Consider the integral 

/ de[i + cos{e - ^^)]e^fc^(i-'=«^(^-^°)) 

Jo 

For large kr the stationary phase points are found from 

^[l-cos(^-^<,)]=sin(^-^„) = 

or ^ = 9o, Oo + n within the range [0, 27r]. Near the first stationary point the integrand 
is dominated by 

When the limits are approximated by (—00, 00), the inegral can be evaluated to give 

J -00 V k,T 

Near the second stationary point the integral vanishes since 1 + cos(6' — 6'o) = 1 — 1 = 0. 
Hence the result (8.5) follows. 



